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APPEAL BRIEF 
L REAL PARTY IN INTEREST 

The real party in interest is the assignee, Telefonaktiebolaget L M Ericsson 
(publ), a Swedish corporation. 

II. RELATED APPEALS AND INTERFERENCES 

A pre-appeal was filed on November 9, 2005, which in accordance with a 
communication from the office on June 15, 2006, did not result in withdrawal of 
the final action. There are no other appeals related to this subject application. 
There are no interferences related to this subject application. 
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III. STATUS OF CLAIMS 

Claims 1, 3-7, 9-28, and 30-44 are pending. Claims 1, 3-7, 9-17, and 20-43 
stand rejected under 35 U.S.C. §103 as being unpatentable over U.S. Patent 
5,317,596 to Ho in view of U.S. Patent 6,597,745 to Dowling. Claims 18, 19, and 
44 stand rejected under 35 U.S.C. §103 as being unpatentable over U.S. Patent 
5,1 17,418 to Chaffee in view of Dowling. Claims 2, 8, and 29 are cancelled. 

IV. STATUS OF AMENDMENTS 
No amendment has been filed after the final action. 

V. SUMMARY OF THE CLAIMED SUBJECT MATTER 

The independent claims 1, 12, 18, 19, 20, and 30 in this case are directed to 
a frequency domain echo canceller, and independent claim 35 recites a method for 
reducing echo at a transceiver. Beneficial example applications of the frequency 
domain echo cancellers and method include multicarrier modulation methods, 
such as DMT (Discrete MultiTone), where digital subscriber lines are employed 
for high speed data transmission and OFDM (Orthogonal Frequency Division 
Multiplexing) based systems. Figure 1 shows a "near-end" DSL modem having a 
transmitter 12 and a receiver 14 connected to the subscriber line's twisted-pair 18 
(telephone line) via a so-called hybrid circuit 16. The transmitted signal is 
transmitted over a communications channel on the subscriber line to a "far-end" 
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DSL modem (not shown), which has its own transmitter and receiver. Because the 
hybrid 16 cannot be exactly matched to the impedance of the subscriber line, the 
transmitted signal "leaks" back into the receiver through the hybrid creating an 
echo signal. An echo canceller in the "near end" station may be used to cancel the 
echo signal. The echo canceller first estimates the echo signal from the 
transmitted signal, and then subtracts the estimated echo signal from the received 
signal. In order to estimate the echo signal, the echo canceller must model the 
echo path channel. (See page 1, line 9-page 2, line 7). 

A DMT signal that is filtered through a physical channel is subjected to 
inter-symbol-interference (ISI) because all practical physical channels have 
"memory," i.e., a non-zero impulse response. This impairment of the physical 
channel causes the transmitted DMT symbols to interfere with each other (ISI). 
Furthermore, the transitions between the DMT symbol cause transients in the 
received signal. These transients cause interference between the carriers in the 
same DMT symbol. This type of interference is known as inter-channel 
interference (ICI), where the "channel" refers to the DMT subchannel. ISI and ICI 
are two distinct phenomena that both cause data distortion. The echo signal 
received at the near-end station, like the data signal transmitted over the channel 
and received at the far-end station, is affected by ISI and ICI, because all practical 
physical echo path channels have a long impulse response (memory). Thus, the 
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echo from one carrier will leak into every other carrier within the same DMT 
symbol and the next transmitted DMT symbol. (See page 3, lines 8-21). 

The frequency domain echo canceller recited in independent claim 20 takes 
into account interchannel interference when estimating an echo signal to be 
removed from a received signal. The frequency domain echo canceller recited in 
independent claim 30 takes into account intersymbol interference and interchannel 
interference when estimating an echo signal to be removed from a received signal. 
For example, see block 50 in Figure 3. Calculating the received echo signal 
completely in the frequency domain is particularly beneficial in a DMT type 
system because the transmitted data is already available in the frequency domain. 
Accordingly, the echo estimate is determined using a frequency domain model of 
an echo path channel that includes the effects of interference like ISI and ICI. 
(See page 4 and page 8, lines 22-page 9, line 2). 

Seven non-limiting example embodiments of a frequency domain echo 
canceller are disclosed and claimed. In a first example embodiment to which 
independent claims 1 and 18 are directed, the frequency domain model of the echo 
path channel is determined using a first matrix of coefficients H, and a second 
matrix of coefficients W,. The first matrix is combined with a current transmitted 
symbol X„ and the second matrix is combined with a previously transmitted 
symbol X/./. The sum of these two combinations is used to estimate the echo 
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signal Y/. The coefficients of the first matrix represent how an echo from a 
currently transmitted frequency domain signal affects a received signal. The 
coefficients of the second matrix represent how an echo from a previously 
transmitted domain symbol affects the received symbol. The coefficients of the 
first and second matrices are adjusted using a difference between the received 
signal and the estimated signal. (See page 9, line 3-page 14, line 20 and Figures 4 
and 5). 

In a second example embodiment, the current transmitted symbol and the 
previously transmitted symbol are partitioned into real and imaginary parts before 
being combined with matrices as described in the first example embodiment. This 
operation reduces the computational complexity by a factor of two compared to 
the first example embodiment. (See page 14, line 21 -page 17, line 9). 

In a third example embodiment to which independent claims 12 and 19 are 
directed, the currently transmitted symbol is combined with a first column vector 
Vj. The previously transmitted symbol is multiplied with a complex exponential 
term, e j2nLklN , in order to compensate for the cyclic prefix, and then subtracted 
from the current transmitted symbol. The resulting signal is combined with a 
matrix Zj . The vector and matrix combinations are summed and used to estimate 
the echo signal to be removed from the received symbol. See equations 40 and 41 . 
This third example embodiment reduces the computational complexity 
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approximately by a factor of two compared to the first example embodiment. (See 
page 17, line 10-page 18, line 17 and Figure 6). Independent method claim 35 is 
directed to both the first and third embodiments. 

In a fourth example embodiment, the operations described in the second 
and third example embodiments are combined resulting in a reduction of 
computational complexity approximately by a factor of four compared to the first 
example embodiment. (See page 18, line 18-page 19, line 6). 

In a fifth example embodiment, when the transmitter of a transceiver has a 
lower sampling rate than its receiver, the received echo signal is interpolated at the 
receiver. The interpolation may also be combined with the operations described in 
the second example embodiment to reduce the computational complexity by a 
factor of four compared to the first example embodiment. (See page 19, line7- 
page20, line 10). 

In a sixth example embodiment, when the transmitter has a higher rate than 
the receiver, the echo signal is decimated at the receiver. The decimation may also 
be combined with the operations described in the second example embodiment to 
reduce the computational complexity by a factor of four compared to the first 
example embodiment. (See page 20, line 1 1 -page 21 line 12). 

In a seventh example embodiment, when the transmitted symbols are not 
aligned in time with the received symbols or frames, an asynchronous echo 
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canceller may be used. The asynchronous echo canceller may be combined with 
previously-described example embodiments. (See page 22, lines 3-21). 

In a DMT type modem, because all of the data is already available in the 
frequency domain, there is no need to do any additional Fourier transform or 
inverse Fourier transform operations, except for the asynchronous echo canceller 
embodiment which requires one extra IDFT. Moreover, some of the coefficients 
in the matrices may be very small and ignored, further reducing the computational 
complexity and memory required. 

VI, GROUNDS OF REJECTION TO BE REVIEWED ON APPEAL 

The two grounds of rejection to be reviewed on appeal are the obviousness 
rejections under 35 U.S.C. §103 based on (1) the combination of Ho in view of 
Dowling, and (2) the combination of Chaffee in view of Dowling. 

VII. ARGUMENT 

A, Ho and Chaffee Lack Multiple Claim Features 

Independent claims 1, 12, 18, 19, and 35 all recite estimating echo in the 
frequency domain using a specific approach. The Examiner admits that Ho and 
Chaffee do not teach estimating the echo signals in the frequency domain using a 
combination of (i) a product of a first matrix (specifically a vector in claim 12) of 
coefficients in the frequency domain and a transmitted symbol and (ii) a product 
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of a second matrix of coefficients in the frequency domain and a previously- 
transmitted symbol. 

More fundamentally, independent claims 1, 12, 18, and 19 recite 
"estimat[ing] in the frequency domain an echo signal" and " removing] in the 
frequency domain the estimated echo signal in the frequency domain from a 
received signal in the frequency domain." Chaffee disclose an echo canceller in 
which the estimated echo is converted to the time-domain and subtracted from the 
received signal. Chaffee does not cancel echoes in the frequency domain. 

Ho relies primarily on echo cancellation in the time domain. See lines 10- 
13 of col. 6: "After the adder 52 subtracts the time-domain portion of the echo, 
e(n), the output is converted by a serial-to-parallel (S/P) converter 54 into a block 
of N real- valued time-domain samples." Then a residual echo is converted from 
the time to the frequency domain so that a residual echo E(f) is removed. Hence, 
the main echo is removed in the time-domain, and then the result is corrected in 
frequency-domain. The reason Ho takes this approach is that Ho does not (and 
can not) remove the echo due to a previous symbol in the frequency domain. 

The Examiner relies on Dowling to remedy the deficiencies in Ho and 
Chaffee. That reliance is misplaced. 
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B. Dowling Is A Pre-Equalizer — Not An Echo Canceller 

Dowling describes a "pre-coder" which is a pre-equalizer — not an echo 
canceller. Dowling precodes a signal before the signal is transmitted in order to 
compensate for transmission channel distortion (not echo) that would otherwise 
distort the signal received at the intended receiver, i.e., the far-end station. The 
hope is that this precoded signal will be received without the far-end station 
having to equalize the signal. In other words, Dowling tries to compensate for the 
channel distortion before transmission over the channel so that the transmitted 
signal is received at the other end of the transmission channel more or less 
undistorted (at least in theory). So Dowling is concerned about the signal received 
at the far-end station. Compensating for the distortion caused by a transmission 
channel is an entirely different problem as compared to canceling an echo. 

In contrast, the claimed echo cancellation is at the near-end station, and the 
echo that the near-end station must deal with is reflected back through an echo 
path to the near-end station. These basic facts would be well understood by those 
skilled in echo cancellation. 1 While Dowling effectively pre-equalizes (precodes) 
the transmitted signal for the transmission channel in order to simplify signal 
processing in the far-end station, the instant claims remove echo in the near-end 



1 For a brief background regarding echo cancellation, please see pages 1 and 2 of the instant specification 
and Figure 1 . 
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station caused by the near-end transmitted signal These are two very different 
operations responding to two different problems. 

C. Compensating For Channel Distortion Does Not Compensate For Echo 

Dowling's precoder does not estimate and remove an echo signal from a 
received signal — let alone do this in the frequency domain. Rather, Dowling 
estimates and "removes" the expected transmission channel distortion from the 
transmitted signal so that distortion introduced by the transmission channel will be 
minimized at the far-end station receiver. Dowling's channel distortion estimate 
does not include echo cancellation. Echo is caused by transmission line 
impedance mismatches, which is different from the channel transfer function. 
Echo is experienced at the near-end station transmitter, and not at the far-end 
station receiver, which is the intended beneficiary of Dowling's precoding. Why 
would Dowling model/estimate an echo that is not part of the channel distortion of 
the signal received at the far-end station receiver? Echo is a problem for the near- 
end station to deal with and not the far-end station. 

The Examiner relies on col. 22, lines 1-3, where Dowling states that 
"communication systems often involve other elements such as echo cancellers 
which may be advantageously merged with the precoder." But here Dowling 
admits that his precoder/pre-equalizer does not perform echo cancellation. If 
echo is to be canceled, Dowling explicitly states that an echo canceller is needed 
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because his precoder does not cancel echo. Thus, the Examiner's premise that 
Dowling's precoder some how removes the echo from the transmitted signal 
cannot stand. 

D. The Ho/Chaffee/Dowling Combinations Fall Short . 

Even if Dowling's precoder were "merged" with Ho's or Chaffee's echo 
canceller, one would not arrive at the claimed echo cancellation. Dowling's 
precoder, by Dowling's own admission, does not perform echo cancellation. 
Hence, in the Examiner's proposed combination, it is just Ho's echo canceller or 
Chaffee's echo canceller doing the echo cancellation, which the Examiner rightly 
admits is not the claimed echo cancellation. So the proposed combination fails. 

To cancel echo, (which is not transmitted over the transmission channel), a 
model of the echo is needed; however, a model of the transmission channel is not. 
On the other hand, Dowling's precoding modifies the transmitted signal using an 
inversion of the transmission channel transfer function (not an echo channel 
transfer function). The hoped for result is the receiver gets a signal that is more or 
less unaffected by the transmission channel transfer function. But only knowledge 
of the transmission channel is needed and not the echo experienced at the near-end 
station. And in the merged combination of Ho and Dowling or Chaffe and 
Dowling, where Dowling's precoder is used along with a separate and distinct 
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echo canceller, the echo is removed from the transmitted signal in a conventional 
way rather than as claimed. 

E. The Ho/Chaffee/Dowling Combinations Would Not Be Made By A 
Person Of Ordinary Skill In The Art 

Notwithstanding the significant deficiencies already noted, assume for the 
sake of argument that there was some provision by Dowling to precode the signal 
to be transmitted for echo cancellation — which Dowling plainly does not do — 
along with precoding for the transmission channel. The echo would be modeled 
and then inverted (call this X) just like the transmission channel is modeled and 
inverted (call this Y). But the consequence is that the inverted echo model X 
negatively affects the inverted channel model Y resulting in distorting the received 
signal even more at the far-end station receiver. For this reason alone, a person 
skilled in this art would not make this modification. 

Moreover, further distorting of the signal at the far-end receiver and the 
need for more sophisticated equalization at the far-end receiver are the very things 
Dowling is trying to avoid. Dowling' s title is "Reduced Complexity Multicarrier 
Precoder," and column 2, lines 41-42 state: "The foregoing indicates a recognized 
but unmet need for a reduced complexity DMT-THP [a DMT-THP is defined in 
column 1, lines 23-24 as "a transmitter-based precoder structure"]." Thus, the 
Examiner requires in the combinations upon which the final rejections are based a 
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modification that renders Dowling inoperable for Dowling's consistently 
articulated and intended purpose. This thwarting of Dowling's primary purpose is 
a clear indicator of an inappropriate obviousness rejection. In re Gordon, 733 
F.2d 900, 902 (Fed. Cir. 1984). 

Consider an alternative combination where Dowling's precoded signal is 
directly applied to the input of Ho ! s echo canceller. In this case, that applied 
signal is still only compensated for transmission channel ISI and ICI. The echo 
ISI and ICI would not be compensated for, but instead would be additionally 
distorted . Accordingly, Ho f s echo canceller cannot be combined with Dowling's 
precoder in the manner proposed by the Examiner to result in the claimed echo 
canceller. Such a proposed combination certainly does not result in estimating a 
echo signal in the frequency domain using the claimed combination of first and 
second matrices with a transmitted symbol and a previously-transmitted symbol. 

These significant obstacles and inconsistencies with trying to combine 
Dowling with Ho and with Chaffee demonstrate that the requisite motivation to 
combine is missing. A proper motivation to combine requires an appreciation of 
the desirability of making the combination. It is not measured by the feasibility of 
making the combination. See Winner Int 7 Royalty Corp. v. Wang, 202 F.3d 1340, 
1349 (Fed. Cir. 2000). In this case, the combination has been shown to be both 
undesirable and infeasible. 

-13- 

1093807 



Appeal Brief 
Lindqvist et al. 
Serial No. 09/584,796 

Still further, the Examiner must show reasons that the skilled artisan, 
confronted with the same problems as the inventor and no knowledge of the 
claimed invention, would select the elements from the cited prior art references for 
the combination in the manner claimed. In re Rouffet, 149 F.3d 1350, 1357 (Fed. 
Cir. 1998). None of Ho, Chaffee, or Dowling reduces the negative affects of ICI 
and ISI in echo cancellation using a frequency domain echo canceller. 
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F. Ho And Dowling Do Not Compensate for Echo ICI or ISI 

The features recited in claims 20 and 30 are not taught in the combination 
of Ho and Dowling. The frequency domain echo canceler in claim 30 "estimate[s] 
the echo in the received signal using a frequency domain model of an echo path 
channel that includes effects of intersymbol interference and inter-carrier 
interference and to subtract the echo estimate from the received signal to provide 
a difference." See also claims 20 and 21. Neither Ho nor Dowling teaches an 
echo canceller that takes into account effects of echo inter-carrier interference 
(ICI) and inter-symbol interference (ISI) in a frequency domain model of an echo 
estimate so that it can be subtracted out from the received echo at the near-end 
station. 

The Examiner assumes that ICI is a form of either ISI or noise. To the 
contrary, as explained in the background of this application, ICI is not the same as 
ISI. Nor is ICI simply "noise." The two articles previously-provided to the 
Examiner and attached in the Evidence Appendix confirm this. 

Another fundamental point ignored by the Examiner is that the ISI and ICI 
in the transmitted signal sent over the transmission channel are not the same as the 
ISI and ICI in the echo . None of the applied references model or compensate for 
the echo ISI and echo ICI. 
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Still further, the rejection of claims 20 and 30 is improper because the 
combination of Ho and Dowling cannot be properly made for the myriad reasons 
already explained. 



Multiple features of the independent claims are not disclosed or suggested 
by the combination of Ho and Dowling or Chaffee and Dowling. Nor is there 
proper motivation to combine their teachings as the Examiner proposes. Each 
missing claim feature and the lack of motivation for each combination is an 
independent ground for reversal. The Board should reverse the outstanding 
rejections. 



VIL CONCLUSION 



Respectfully submitted, 



NIXON & VANDERHYE P.C. 




JRL/maa 
Enclosures 



Reg. No. 33,149 
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IX, CLAIMS APPENDIX 

1 . An echo canceller for use in a transceiver, comprising: 
first electronic circuitry configured to estimate in the frequency domain an echo 
signal, and 

second electronic circuitry configured to remove in the frequency domain the 
estimated echo signal in the frequency domain from a received signal in the frequency 
domain, 

wherein the first electronic circuitry is further configured to estimate the echo 
signals in the frequency domain using a combination of (i) a product of a first matrix of 
coefficients in the frequency domain and a transmitted symbol and (ii) a product of a 
second matrix of coefficients in the frequency domain and a previously-transmitted 
symbol. 

3. The echo canceller in claim 1, wherein transmitted signals corresponding to the 
transmitted symbol and the previously-transmitted symbol are real- valued, and wherein 
the transmitted symbol and the previously-transmitted symbol are divided into real and 
imaginary parts before being combined with the respective matrix to reduce 
computational complexity. 

4. The echo canceller in claim 1, wherein the coefficients of the first matrix 
represent how an echo from a currently transmitted frequency domain signal affects a 
received signal. 
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5. The echo canceller in claim 4, wherein the coefficients of the second matrix 
represent how an echo from a previously transmitted frequency domain signal affects the 
received signal. 

6. The echo canceller in claim 1, wherein the first electronic circuitry is further 
configured to adapt the coefficients of the first matrix and the second matrix using a 
difference between the received signal and the estimated echo signal. 

7. The echo canceller in claim 6, wherein the first electronic circuitry is further 
configured to adapt the coefficients using a least mean squares algorithm. 

9. The echo canceller in claim 1, wherein the transceiver is a discrete multitone 
(DMT) transceiver. 

10. The echo canceller in claim 1, wherein the first and second matrices are Nx 
N matrices, where N is a number of symbol samples. 

1 1 . The echo canceller in claim 1 , wherein a vector corresponding to a transmitted 
frequency domain symbol, a vector corresponding to a received frequency domain signal, 
and a vector corresponding to an estimate of the echo symbol are all Hermitian 
symmetric. 

12. The echo canceller for use in a transceiver, comprising: 

first electronic circuitry configured to estimate in the frequency domain an echo 
signal, and 
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second electronic circuitry configured to remove in the frequency domain the 
estimated echo signal in the frequency domain from a received signal in the frequency 
domain, 

wherein the first electronic circuitry is further configured to estimate the echo 
signals in the frequency domain using a combination of a product of (i) a vector of 
coefficients in the frequency domain and a transmitted symbol and (ii) a product of a 
matrix of coefficients in the frequency domain and a compensated, previously- 
transmitted symbol . 

13. The echo canceller in claim 12, wherein the first electronic circuitry is further 
configured to divide the transmitted symbol and the previously-transmitted symbol into 
real and imaginary parts before combining them respectively with the vector and the 
matrix to reduce computational complexity. 

14. The echo canceller in claim 12, wherein a compensation factor used to 
compensate the previously-transmitted signal is a complex exponential term. 

15. The echo canceller in claim 14, wherein the transceiver is a discrete multitone 
(DMT) transceiver and the compensation factor compensates for a cyclic prefix 
associated with the previously-transmitted signal. 

16. The echo canceller in claim 1, wherein when a transmitter of the transceiver 
has a lower sampling rate than a receiver of the transceiver, the echo signal is interpolated 
at the receiver. 
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17. The echo canceller in claim 1, wherein when a transmitter of the transceiver 
has a higher sampling rate than a receiver of the transceiver, the echo signal is decimated 
at the receiver. 

18. An echo canceller for use in an asynchronous transceiver configured to cancel 
an echo signal, comprising: 

a first matrix of coefficients; 

a second matrix of coefficients; and 

electronic circuitry configured to use a combination of (i) a product of the first 
matrix and a currently-transmitted symbol and (ii) a product of the second matrix and a 
previously-transmitted symbol to estimate an echo signal in the frequency domain, to 
transform the estimate of the echo signal into the time domain, and to remove the 
transformed estimate from a received signal in the time domain. 

19. An echo canceller for use in an asynchronous transceiver configured to cancel 
an echo signal, comprising: 

a vector of coefficients in the frequency domain; 
a matrix of coefficients in the frequency domain, 

electronic circuitry configured to use a combination of (i) a product of the vector 
and a currently-transmitted symbol and (ii) a product of the matrix and a compensated, 
previously-transmitted symbol to estimate an echo signal in the frequency domain, to 
transform the estimated echo signal into the time domain, and to remove from a received 
signal in the time domain. 
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20. An echo canceller for use in a transceiver canceling an echo from a received 
signal in the frequency domain including circuitry configured to determine an estimate of 
the echo in the received signal using a frequency domain model of an echo path channel 
that takes into account effects of inter-carrier interference and to subtract the echo 
estimate from the received signal. 

21 . The echo canceller in claim 20, wherein the echo canceller is used in a 
discrete multitone (DMT) type transceiver and the frequency domain model takes into 
account intersymbol interference . 

22. The echo canceller in claim 20, wherein the frequency domain model includes 
a first set of values that models how an echo from a currently transmitted frequency 
domain symbol distorts the received signal and a second set of values that models how an 
echo from a previously transmitted frequency domain symbol distorts the received signal. 

23. The echo canceller in claim 22, wherein the first set of values is a first 
complex matrix and the second set of values is a second complex matrix. 

24. The echo canceller in claim 22, wherein the first set of values is a column 
vector and the second set of values is a matrix. 

25. The echo canceller in claim 24, wherein the matrix is combined with a 
difference between the currently transmitted symbol and a product of the previously 
transmitted symbol and a compensating factor. 
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26. The echo canceller in claim 22, wherein the transmitted symbol and the 
previously transmitted symbol are divided into real and imaginary parts before being 
combined with the first and second sets of values, respectively. 

27. The echo canceller in claim 20, wherein when a transmitter of the transceiver 
has a lower sampling rate than a receiver of the transceiver, the echo signal is interpolated 
at the receiver. 

28. The echo canceller in claim 20, wherein when a transmitter of the transceiver 
has a higher sampling rate than a receiver of the transceiver, the echo signal is decimated 
at the receiver. 

30. A frequency domain echo canceller for use in a transceiver canceling an echo 
from a received signal in the frequency domain including circuitry configured to 
determine an estimate of the echo in the received signal using a frequency domain model 
of an echo path channel that includes effects of intersymbol interference and inter-carrier 
interference and to subtract the echo estimate from the received signal to provide a 
difference. 

3 1 . The echo canceller in claim 30, wherein the frequency domain model includes 
a first set of values that models completely in the frequency domain how an echo from a 
currently transmitted frequency domain symbol distorts the received signal and a second 
set of values that models completely in the frequency domain how an echo from a 
previously transmitted frequency domain symbol distorts the received signal. 
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32. The echo canceller in claim 3 1, wherein transmitted signals corresponding to 
the currently and previously transmitted frequency domain symbols are real-valued. 

33. The echo canceller in claim 32, wherein the currently transmitted symbol, the 
previously transmitted symbol, the received signal, and the difference are vectors having 
Hermitian symmetry. 

34. The echo canceller in claim 3 1 , wherein the difference is used to adjust the 
first and second set of values. 

35. A method for reducing an echo at a transceiver comprising: 

(a) combining in the frequency domain a currently-transmitted symbol with a 
first vector or matrix of coefficients in the frequency domain resulting in a first 
combination; 

(b) combining in the frequency domain a previously-transmitted symbol with a 
second matrix of coefficients in the frequency domain resulting in a second combination; 

(c) combining the first and second combinations in the frequency domain to 
estimate the echo in the frequency domain; and 

(d) using the estimated echo to reduce the echo in a signal received at the 
transceiver. 

36. The method in claim 35, further comprising: 

determining a difference between the received signal and the estimated echo, and 
adjusting the first and second set of values using the difference. 
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37. The method in claim 35, wherein the first vector or matrix corresponds to a 
first matrix of coefficients . 

38. The method in claim 35, wherein the first vector or matrix corresponds to a 

column vector of coefficients . 

39. The method in claim 38, wherein the combining step (b) includes: 
multiplying the previously transmitted symbol by a compensation factor to 

produce a product; 

subtracting the product from the currently transmitted symbol; and 
combining a result of the subtracting with the matrix. 

40. The method in claim 35, wherein when a transmitter of the transceiver has a 
lower sampling rate than a receiver of the transceiver, the method further comprising: 

interpolating the echo signal. 

41 . The method in claim 35, wherein when a transmitter of the transceiver has a 
higher sampling rate than a receiver of the transceiver, the method further comprising: 

decimating the echo signal. 

42. The method in claim 35, wherein a vector corresponding to a transmitted 
frequency domain symbol, a vector corresponding to a received frequency domain signal, 
a vector corresponding to an estimate of the echo symbol are Hermitian symmetric. 

43. The method in claim 35, further comprising: 
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dividing the currently transmitted symbol and the previously-transmitted symbol 
into real and imaginary parts before the combining steps (a) and (b), respectively, to 
reduce computational complexity. 

44. The method according to claim 35, wherein the transceiver is an asynchronous 
transceiver, the method further comprising: 

transforming the estimated echo into the time domain, and 

removing in the time domain, the estimated echo signal from the received signal 
on a sample-by-sample basis. 
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X. EVIDENCE APPENDIX 

Two articles are attached that were submitted to the Examiner on June 15, 
2005 and are referred to in the argument section above: "A Novel Intercarrier 
Interference Cancellation Approach in OFDM based on BSS," by Liu et al and "Residual 
Frequency Offset Correction for Coherently Modulated OFDM Systems in Wireless 
Communication," by Abhayawardhana et al. 

XI. RELATED PROCEEDINGS APPENDIX 

There is no related proceedings appendix. 
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Abstract — Orthogonal Frequency Division Multiplex ( OFDM) systems are 
very sensitive to frequency offset caused by tuning oscillator instabilities and 
Doppler shifts induced by the channel. The Schmidl and Cox Algorithm 
(SCA) is quite robust in estimating the frequency offset for systems with large 
OFDM symbol lengths. It uses two OFDM symbols for training with the first 
one having two identical halves. The frequency offset is estimated by correlat- 
ing a received sequence of samples equal to half the OFDM symbol length with 
the following received samples. The effect of Additive White Gaussian Noise 
(AWGN) in the estimation process is mitigated only if the number of samples 
used in the correlation, and hence the OFDM symbol size is large. However to 
be successfully applied to Broadband Fixed Wireless Access (BFWA) systems, 
OFDM should perform well even with smaller symbol lengths. In this paper 
we present the Residual Frequency Offset Correction Algorithm (RFOCA), 
which uses the SCA for initial frequency offset acquisition and follows it with 
a stage which reduces the initial residual frequency offset by tracking the 
phase of the decoded data. We show through simulation that the RFOCA 
achieves an error variance that is many orders of magnitude lower than at 
the end of the acquisition stage using the SCA alone. 



I. Introduction 

Various solutions have been suggested as contenders to over- 
come the challenge of transmitting data to the subscriber at high 
data rates demanded by future broadband applications, for ex- 
ample Broadband Fixed Wireless Access (BFWA), Digital Sub- 
scriber Lines (DSL), Cable Modems and Satellite Communica- 
tion. This paper will address BFWA, owing to its advantages 
of ease of deployment and cost Orthogonal Frequency Divi- 
sion Multiplex (OFDM) has been the most popular physical layer 
standard for BFWA systems due to its robustness and efficiency 
mainly because it uses a simple Fast Fourier Transform (FFT) 
of finite length N for modulation and demodulation. Among 
the standards that have placed their confidence in OFDM are 
HIPERLAN-2 [1] and IEEE 802.1 1 -a [2]. 

The orthogonality of the consecutive OFDM symbols is main- 
tained by appending a length v cyclic prefix (CP) at the start of 
each symbol [3]. The CP is obtained by taking the last v samples 
of each symbol and consequently the total length of the transmit- 
ted OFDM symbols is (N + v) samples. For each OFDM symbol 
to be independent and to avoid any Inter Symbol Interference (ISI) 
or Inter Carrier Interference (ICI), the length of the Channel Im- 
pulse response (CIR) should be less than v + 1 samples. Hence the 
distortion caused by the CIR only affects the samples within the 
CR The receiver discards the CP and takes only the last N samples 
of each OFDM symbol for demodulation by the receiver FFT. In 
such a case, the linear convolution of the transmitted signal with 
the CIR is converted into a circular one. Consequently, for coher- 
ently modulated OFDM systems the effects of the CIR can then 
be equalized by an array of one-tap Frequency Domain Equalizers 
(FEQ) following the FFT. This is because the frequency selective 
fading channel can be approximated by a sum of flat fading chan- 
nels, provided the number of subchannels is large enough. Some- 
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Fig. 1. Block Diagram of the OFDM System 



times a Time Domain Equalizer (TEQ) may be used to shorten 
the effect of the CIR [4]. Figure 1 shows the block diagram of 
the system, where P/S and S/P represent the Parallel to Serial and 
Serial to Parallel functions respectively. 

Unfortunately, OFDM has been proven to be very sensitive to 
frequency offset and phase noise caused by tuning oscillator in- 
accuracies or Doppler shifts induced by the channel. The usual 
deployment scenario for a BFWA system employs a point to mul- 
tipoint configuration. Here, a single base station (BS) transmits 
data in short bursts to many subscriber units (SUs) placed at the 
user locations. The economic feasibility of a BFWA deployment 
depends heavily on the cost of the SUs. The use of low cost com- 
ponents, particularly the oscillators, is a major issue since their 
accuracy and stability are directly related to cost. Hence it is im- 
perative that if OFDM is to be used for BFWA transmission, it 
should be able to operate effectively using oscillators with only 
moderate performance and cost. 

An OFDM system effectively consists of N sinusoidal subcar- 
riers with frequency spacing 1/T, where T is the active sym- 
bol period of each subcarrier. The fcth subcarrier will thus be at 
fk = fo + k/T, where /o is a reference frequency. Without loss 
of generality, we can assume fo = 0. The modulated subcarri- 
ers overlap spectrally, but since they are orthogonal over a symbol 
duration, they can be easily recovered as long as the channel does 
not destroy the orthogonality. An unwindowed OFDM system has 
rectangular symbol shapes. Hence, in the frequency domain the 
individual subchannels will have the form of sine functions where 
the first sidelobe is only some 13 dB below the main lobe of the 
subcarrier. In the event of a frequency offset between the trans- 
mitter and the receiver the subchannels will not align. The effect 
of frequency offset is two fold, it will reduce the signal power of 
the desired subchannel outputs and also introduce ICI further in- 
creasing the interference levels. The ICI power can be significant 
even for small frequency offsets in OFDM due to the high sidelobe 
power of the subchannels. Thus the reduction of frequency offset 
is critical in OFDM systems. This paper is organised as follows, 



section II analyses the effect of frequency offset in OFDM analyt- 
ically. Section III briefly introduces the Schmidl and Cox Algo- 
rithm (SCA) [5] and section IV introduces the novel Residual Fre- 
quency Offset Correction Algorithm (RFOCA) that compliments 
the SCA and improves the performance. Section V shows results 
obtained by computer simulation and section VI concludes and 
proposes future work. 

II. Effect of Frequency Offset in OFDM 

All analysis and simulations in this paper are performed in the 
digital complex baseband domain. The nth sample of the mth 
OFDM symbol generated by the Inverse FFT (IFFT) at the trans- 
mitter is 



N ^ 



A m , k e j2 *ft,Q<n<N-l 



(1) 



A mi k is the data symbol modulated on to the fcth sub- 
carrier of the mth OFDM symbol. The data is con- 
verted into a serial sequence, then the CP of length v is 
added. Thus the mth transmitted OFDM symbol is s(m) = 
[5 m ,iv- v ,..,5 m ,jv-i J 5 m) o,..,Sm t iV_i] r . The receiver discards 
the samples of the CP and uses the samples of the active OFDM 
symbol for decoding. We assume a finite length CIR with N h 
samples, h = [ho, .., ^-i] T , where v > N h - 1. We assume 
that N is large enough such that the frequency selective channel 
is divided into contiguous flat fading subchannels. If we assume 
that correct frame and timing synchronisation is achieved, then 
the received sequence after stripping the CP can be expressed as, 



1 V Jfe=0 



(2) 



for 0 < n < N - 1. Here H k is the transfer function of the 
channel at the subchannel index k, w n is the Additive White Gaus- 
sian Noise (AWGN) term and e is the frequency offset relative to 
the intercarrier spacing. The symbol after FFT demodulation is 



n=0 



(3) 



Substituting (2) into (3) allows the output of the Ith subchannel 
to be expressed as the sum of three components. 



Y m ,i = (A roi ,ff|) 



sin(7re) 
iVsin(7re/iV) 



j™-*- + jj + W t (4) 



The first term on rhs of equation 4, denoted Ri hereafter, shows 
that the desired term A m ,iHi experiences amplitude reduction and 
phase distortion. As N » 7re, N sin(7re/iV) can be approximated 
by 7re. Hence the degradation of the required term vanishes when 
the received data does not experience any frequency offset. In 
practice, the effect of the channel transfer function Hi is removed 
by the FEQ following the FFT. Note that the FEQ is trained by 
transmitting a known training symbol occupying the OFDM sym- 
bol number N tt A = [A Ntt0 , .., A NttN -i) T . Usually the training 
symbol is placed at the beginning of the frame. The ratio of the 
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Fig. 2. Analysis of CTIR vs Relative Frequency Offset, e 

decoded samples of the training symbol 1jv (j j, to the locally gen- 
erated copy of Aw t where 0 < / < N - 1, is used to calculate 
the coefficients for the FEQ. 

The second term in equation (4) is the ICI and the third term is 
due to the AWGN. The ICI term is given as, 



JV-l 



I, = y £ {Am;kHk) Nsm«(k + e- 

hm,0 X 



l)/N 



(5) 



If we assume that the transmitted data is zero mean and uncor- 
rected (i.e. E[A m ,k] = 0 and E[A mJ A* mJt ] = \A\ 2 6 ltk ), then 
E[I{\ = 0. Hence the ICI power is, 



N-l 



£[|/H 2 ] = M| 2 [sin7re] 2 £ 
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(6) 
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To appreciate the severity of frequency offset on the Carrier 
to Interference Ratio (CTIR) in an OFDM system, we define 
CTIR = #[|^| 2 ]/£[|^| 2 ]- Assuming the channel transfer func- 
tion, H is unity and with zero AWGN, (i.e. Wi = 0), then 



CTIR = 



Z^-o 1 l/[7Vsin7r(fc + e- 
[ATsin(7re/A0] 2 



O/Af 



(7) 



A similar analysis on the effect of frequency offset on OFDM 
can be found in [6] and [7]. Figure 2 shows CTIR in decibels 
as a function of e. Simulated values for N = 64 and 256 using 
QPSK data mapping are also shown. It shows that the CTIR drops 
significantly when relative frequency offset, e increases beyond 
0.1. An important observation is that the CTIR is independent of 
N. 

III. Schmidl and Cox Algorithm (SCA) 

One of the more robust schemes to estimate both frame syn- 
chronisation and frequency offset is the SCA. It uses two training 
symbols with the first one having a repetition within half a sym- 
bol period. Frame synchronisation is achieved by searching for a 
training symbol with two identical halves. If L = N/2, the sum 



of L consecutive correlations between pairs of samples spaced L 
apart is found as, 

n=0 

The output of the correlator, as given by (8) will reach a peak 
at the end of the CP of the first training symbol. If v > Nk - 1 
and p =s v— {N h - 1), the peak is actually maintained for p 
samples just before the end of the CP of the first training symbol. 
Thus the correlator output will take the form of a plateau. This is 
because the last few samples of the CP, equal to the length of p, 
are not corrupted by the CIR. Timing synchronisation is achieved 
by locating the end of this plateau, denoted by d opt . Since the 
CP ensures that there is no Inter Block Interference (IBI) between 
OFDM symbols, the effect of the CIR in the first training symbol 
is cancelled if the conjugate of a sample from the first half is mul- 
tiplied by the corresponding sample from the second half. Hence 
the phase difference between the two halves of the first training 
symbol is caused by the frequency offset, <j> = 7re. It can be esti- 
mated as^ = l/pEn=o£ p ( d °pt- n )' Ife < 1 there is no phase 
ambiguity in 4> and the frequency offset can be estimated as 

csca - <j>fa ( 9 ) 

In order to resolve potential ambiguity use can be made of the 
second training symbol, as detailed in [5]. The effect of AWGN 
on the partial estimation given by equation (9) is insignificant only 
if L, and hence the number of FFT points JV, is large. The SCA 
performs well for OFDM systems with JV in excess of 1000 [5]. 
For BFWA systems, data is transmitted in short bursts, particu- 
larly in the uplink. In this situation, it would be wasteful to use 
an OFDM system with large JV. Besides large values of JV will 
give rise to the additional problems of high peak-to-average power 
ratio and will also introduce latency which reduces protocol effi- 
ciency. Typically BFWA systems utilise lower values of JV such 
as 64 to 256 as presented in [1] and [2]. When the SCA is used for 
systems with lower values of JV, the estimate i is not sufficiently 
accurate. This error results in a residual frequency offset that ro- 
tates the received constellation at a reduced rate, but one that is 
still significant enough to cause bit errors in coherently demodu- 
lated OFDM systems. In this paper we present the Residual Fre- 
quency Offset Correction Algorithm (RFOCA) that compliments 
the SCA by following it with a tracking function that continuously 
compensates for the residual frequency offset error. 

Other approaches that are used for frequency offset correction 
include for example the use of null symbols [8]. The use of a null 
symbol is not practical in a burst mode system as the transmission 
will always be a null while no data is transmitted. Other sugges- 
tions utilise the correlation between the cyclic prefix and the last 
few samples of the OFDM symbol [9] and the seif-cancellation 
schemes presented in [10] and [7], where the same data is sent in 
more than one subchannel. Both schemes are not very practical 
when JV is a low value with the latter two approaches having the 
additional problem of reducing the number of useful subchannels 
by at least half. Note that the RFOCA does not utilise any pi- 
lots nor null symbols yet will be shown to effectively eliminate 
the residual offset. The proposal in [1 1] also attempts to track the 
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Fig. 3. Residual Frequency Offset Correction Algorithm (RFOCA) 

residual frequency offset. Here, the residual offset correction fac- 
tor is updated at the end of each OFDM sy mbo 1 . However our tests 
indicate that the estimate is not very accurate when only one sym- 
bol is used, hence we update the residual offset correction factor 
only after analysing a block of N w , (N w > 1) OFDM symbols. 

IV. Residual Frequency Offset Correction 
Algorithm (RFOCA) 

In the proposed approach, an initial frequency offset acquisition 
is made using the SCA. We propose to estimate the residual fre- 
quency offset, e by tracking the rate of phase change at the FEQ 
outputs, Y m j as shown in Figure 3. However, the estimate of e 
could be seriously affected by subchannels with a low SNR result- 
ing from spectral nulls in the channel response Hi, An estimate 
of Hi is made by taking the ratio of the transmitted and decoded 
output of the second training symbol of the SCA. Hence we only 
select those subchannels with \H t | above a certain threshold. We 
call this subset of subchannels d C [0, ..JV - 1]. The criteria ap- 
plied is to select subchannels with \Hi\ in excess of a standard 
deviation above the mean. For symbol m, the outputs of these 
subchannels Y m ^ are sent through a slicer to obtain F m ,d, where 

ded. 

For symbol m the phase errors between Y m j and Y m ,d, namely 
0m,rf are found as, 

0 m4 = ZY m , d - LY mA (10) 

The phase errors are stored in a buffer, b for a block of N w OFDM 
symbols. Note that only the selected subchannels will be stored. 
Hence if the number of subchannels selected in d is N d then the 
length required for the buffer is only N d N w . The Maximum Like- 
lihood (ML) estimate of the residual frequency offset e is calcu- 
lated as the gradient of the values stored in the buffer. Note, in 
this paper we have assumed QPSK data mapping though the tech- 
nique could be extended to other constellations. Initially, i can be 
quite high and if the values of the phase errors, 0 m , rf > 7r/4, then 
it results in decoding errors being produced by the slicer which 
subsequently results in phase wrapping of 0 m)rf at ±t/4. Hence 
it is imperative that we first unwrap 6 m ^ before the calculation of 
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Fig.4. Phase Unwrapping Algorithm 

I. Although accurate phase unwrapping algorithms are available, 
they are very complex. Since we are only interested in the phase 
gradient and not the exact phase values, a simpler scheme was 
selected for the unwrapping of the phase [12]. For purposes of 
clarity, we denote the ith sample of the wrapped phase and the un- 
wrapped phase as $i and 0 it respectively. The phase unwrapping 
algorithm can be expressed as 

Si = + aSAW(0i - Oi-i) (11) 

where SAW(.) is a sawtooth function that limits the output to 
±7r/4 and a is a parameter that controls the variance of the un- 
wrapped phase. The algorithm is illustrated in Figure 4. 

The error due to frequency offset is compensated by multiply- 
ing the pre-FFT received symbols, r m}Tl with e jnim/N y where 
in general e m is the frequency offset correction factor for sym- 
bol m, which is initialised to i S cA at the acquisition stage. At 
symbol number qN wt e q N„ is updated according to i qNw ~ 
e q N w -i + £ q N v , where q is an positive integer. For OFDM" sym- 
bols" [fa - l)N w ,qN w - 1], € (ff -i)Ar w is used as the frequency 
offset correction factor. In other words i m is updated once every 
N w OFDM symbol blocks. 

The choice of N w is critical as a large value will cause the up- 
dating of e m with too high a latency. However after a few updates 
the values of e are quite small, hence to get a better estimate in the 
presence of AWGN a larger value of N w will reduce the variance. 
The RFOCA can be summarised as follows. 

1. Do an initial frequency offset acquisition iscA using the SCA. 

2. Select the subset of subchannels with magnitudes exceeding 
the defined threshold in the channel transfer function, d, at the 
start of the burst. 

3. Obtain Y m , d from the FEQ output and then Y m , d by use of a 
sheer for each symbol m, where d E d. Calculate and store d m ^ 
for a block of N w OFDM symbols. 

4. Find the unwrapped phase 8 from the wrapped phase 9. Find 
the gradient of namely e. This process is done once every N w 
symbols. 

5. Once every N w OFDM symbol block, calculate the new fre- 
quency offset correction factor i q N w . 

V. Simulation Parameters and Results 

OFDM systems with N = 64, 128, 256 have been simulated at 
a sampling rate of 20 MHz with a guard interval equal to 20 sam- 
ples, thus the subcarrier spacings are approximately 312 kHz, 156 
kHz and 78 kHz respectively. QPSK mapping for all subchannels 
has been employed and all the subchannels are used. A burst of 
320000 data bits is assumed to be transmitted. Each data point 
in the simulation results is obtained by averaging over 500 such 
bursts. In order to test the performance of the RFOCA with fre- 
quency offset alone, we have assumed perfect symbol and frame 
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Fig. 5. Comparison of Error Variances of SCA and RFOCA in AWGN for N 
64, 128,256 and e = 0.5 



synchronisation for the simulated coherent OFDM systems. A 
more robust symbol synchronisation algorithm that virtually guar- 
antees perfect frame synchronisation will be the subject of a fu- 
ture publication. After much testing, the values selected for N w 
were 250, 250, 125 for N = 64, 128, 256, respectively. Figure 5 
shows the comparison of error variances at the end of the acquisi- 
tion stage using the SCA and the tracking stage using the RFOCA 
for A T = 64, 128, 256 as a function of the Signal to Noise Ratio 
(SNR), with e = 0.5. It shows that the RFOCA error variance 
is many orders of magnitude lower than at the end of the acqui- 
sition stage. There appears to be a threshold effect in RFOCA 
at around an SNR of 16 dB, below which the error variance in- 
creases rapidly. It is mainly due to the errors present in the se- 
lected subchannels of d, which subsequently create estimation er- 
rors in the residual offset, i m . Figure 6 shows the performance 
of the RFOCA in terms of Bit Error Rate (BER) vs SNR again in 
AWGN. The BER is virtually zero at values of SNR above 16 dB 
but suddenly increases below L6 dB. This result is a direct con- 
sequence of the threshold effect observed with the previous error 
variance results. It is observed that at an SNR of 16 dB only a few 
out of the total of 500 received bursts give rise to errors. How- 
ever, these errors are significant giving rise to a high overall BER. 
Such occurances became more common when the SNR is reduced 
below 15 dB, thus increasing the BER very rapidly. 

Appropriate models for BFWA channels are in the process of 
being defined. The Stanford University Interim (SUI) channels 
comprise 6 models for 3 different terrain conditions [13]. All 
of them are simulated using 3 taps, each having either Rayleigh 
or Ricean amplitude distributions. The channel is assumed to be 
wide-sense stationary uncorrelated scattering (WSSUS) and each 
tap of the CIR is modeled as hi - fte J '*\ where the amplitude Si 
and the phase are selected independently [14]. The burst takes 
less than 10 ms to transmit at the selected sampling frequency, 
consequently the channel is assumed constant for the duration of 
each burst. We have selected the SUI-2 channel model, pertaining 
to terrains with low tree densities and with antennas having a di- 
rectivity of 30 degrees at the SU and 120 degrees at the BS. The 
channel is characterised by a RMS delay spread of 0.2 fis. 

Figure 7 allows a comparison of the error variances of the SCA 
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Fig. 6. Performance of RFOCA in AWGN for N = 64, 128, 256 and e = 0.5 
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Fig. 7. Comparison of Error Variances of SCA and RFOCA with AWGN and 
SUI-2 CIR for N - 64, 128, 256 and e = 0.5 

and the RFOCA for an OFDM system with e = 0.5 in the SUI- 
2 channel. A new channel in accordance with the SUI-2 profile 
is randomly generated for the transmission of each burst. The 
RFOCA error variance is seen to be many orders of magnitude 
lower than at the end of the acquisition stage despite the presence 
of the SUI-2 channel. 

Figure 8 shows the performance of the RFOCA in terms of BER 
vs SNR for the SUI-2 channel and e = 0.5. The threshold effect 
now occurs at an SNR in the region of 20 dB owing to the effect 
of the SUI-2 channel. 

VI. Conclusion 

We have presented a Residual frequency Offset Correction Al- 
gorithm (RFOCA), that compliments the frequency acquisition 
process performed by the Schmidl and Cox Algorithm (SCA). The 
RFOCA continuously tracks and compensates for the residual fre- 
quency offset that is present after the acquisition stage. We have 
applied the RFOCA for AWGN and BFWA channels, namely the 
SUI-2 channel profile, and have shown a significant reduction in 
the error variance brought about by the use of the RFOCA fol- 
lowing the initial acquisition stage. Although the algorithm may 
appear to suffer from a threshold effect it still gives a significant 
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Fig. 8. Performance of RFOCA with AWGN and SUI-2 CIR for N = 
64,128,256 and e = 0.5 

performance advantage at realistic signal to noise ratios. In the 
future we hope to address the threshold effect by improving the 
phase unwrapping and investigating other criteria for the selection 
of the subchannels. We also hope to improve the performance of 
the SCA in a frame synchronisation role. 
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Abstract- Orthogonal Frequency Division Multiplexing (OFDM) is widely applied in wireless 
communication systems nowadays. In practice, there are frequency differences between the local 
oscillators in the transmitter and in the receiver. Due to these frequency offsets, the sub-carriers are 
not orthogonal. This leads to InterCarrier. Interference (ICI) between the sub-carriers and degrades 
the system performance severely. 

In this contribution, we propose a new ICI cancellation method based on Blind Source Separation 
(BSS). The relative gradient algorithm is employed to produce the separating matrix. Also, a 
technique to solve the permutation ambiguity is developed. The algorithm's interference cancellation 
is maintained for a wide range of interference conditions. Computer simulations are given, which 
confirm the effectiveness of the proposed compensation technique. 

1. Introduction 

In OFDM communication systems, a wideband source signal is partitioned into a number of the 
narrow sub-signals which are transmitted simultaneously employing the orthogonal sub-carriers. 
OFDM has various applications in wireless communications, such as digital audio, digital TV [1], and 
broadband satellite communication [2]. One major limitation of OFDM in many applications is its 
sensitivity to frequency differences between the local oscillators in the transmitter and in the receiver. 
These frequency offsets cause crosstalk between the sub-carriers, namely, ICI. ICI increases the 
Symbol Error Rate (SER) of the demodulated signal and degrades the performance of the receiver 
[3,4]. 

Zhao and Haggman have proposed ICI self-cancellation scheme to mitigate the ICI effects [5]. 
The resulting system is less bandwidth efficient than the normal OFDM scheme [5,6]. 

In this paper, an InterCarrier Interference Cancellation approach based on BSS (ICI/BSS) is 
proposed. In addition, a technique is given to solve the permutation ambiguity inherent in BSS type 
algorithms, which makes ICI/BSS practical in practice. Another advantage of ICI/BSS is that no 
prior information is needed. Also, no training sequences are required, which leads to efficient 
utilization of the bandwidth. 

This paper is organized as follows. For a perfect Nyquist channel, the model of ICI resulting 
from the frequency offsets is described in detail in Section 2. ICI/BSS employing the relative 
gradient algorithm is given in Section 3. The technique to solve the permutation ambiguity is also 
given in Section 3. In Section 4, simulation results are given. Conclusions are drawn in Section 5. 

2. Mathematical Formulation 

The complete block diagram of the transceiver for OFDM system is given [I]. In this paper, the 
analysis considers only the impairments due to the frequency offsets. Thus, the forward error correct 
coding bock is not included. The frequency offsets alone do not cause InterSymbol Interference (ISI). 



Usually, a cyclic prefix is used to eliminate the ISI in OFDM. The use of a cyclic prefix is not 
considered in this analysis too. 

The structure of the simplified OFDM system, used to derive the ICI mathematical model, is 
shown in Fig. 1. In OFDM, interleaving is applied to correct bit errors, which occurs in burst rather 
being randomly scattered. A common used interleaving scheme is block interleaving, where input 
bits are written into a block column by column and read out row by row. The length of the block's 
column equals several times of the symbol length [7]. After interleaving, the serial source signal s(n) 
is independent at the time instant n. In general, s(n) is assumed to be complex-valued, zero-mean, 
stationary, and nongaussian distributed. The symbol time of s(n) is T, In Fig. 1 , the parallel source 
signals Sj(m) (i=0,l,...,N-l) are given by 

s i (m) = s(iN + m) ( l > 

where N is the number of the sub-carriers. 

The symbol time of Si(m) is T p , which is related to T s as 

T„ = NT, W 

Since s(n) is white in the time domain, the signals Si (m) are considered to be mutually 
independent. 

Then, N points Inverse Discrete Fourier Transform (IDFT) is applied to modulate Sj (m) into the 
orthogonal sub-carriers. The modulated orthogonal signals b k {m) (k=0,l,...,N-l) are given by 

where k is the IDFT index. . . 

The signal b(n) is obtained from b k (m) by parallel to serial operation. The transmitted signal 

s-r(t) is given by 

s T (/) = exp(y 2nf c t ) £ b(n)g T (t-nT,) ( 4 ) 

where f c is the local oscillator frequency, and gT (t) is the impulse response of the low-pass filter, in the 
transmitter. 

To simplify the following analysis, the additive white gaussain noise n(t) is ignored m the 

derivation. The received signal s R (t) is given as 

s R (t) = s T (t)®h(t) 0>) 

where h(t) is the impulse response of the channel and ® denotes convolution. 

In the receiver side, the local oscillator frequency is / c + 4f . which is different from the 
transmitter side by A/. The phase delay between the local oscillators in the transmitter and in the 
receiver is ignored. Thus, after the demodulation by / c + A/and the low-pass filter g R (t), in the 
receiver, the signal c(/) is given by 

c(t ) = exp(-;2«V"r ) £ b(n)g T (r - nT) ® MO ® g» (0 

(6) 

= exp(-;2/rA/?)£ b(n)p(t-nT) 

where g R (t) is the impulse response of the low-pass filter in the receiver and p(t) is the combined 
impulse response of the transmitter filter, the channel filter, and the receiver filter. 

Assuming that p(t) satisfies the Nyquist pulse-shaping criterion for samples taken at the optimum 
instants of intervals T s [6]. After ADC and serial to parallel block in the receiver, the signals are 



given by 

c k (m) = exp{j2x*fl:T,)b k (m) - cxp(^)b k (m) (7) 
where the frequency error is defined 

. e - NAfT s = AJT p (8) 

N points Discrete Fourier Transform (DFT) is applied to the signals, c k (m), to produce the 
demodulated signals 

*,(") = X c *< m ) ex K^^) 



**0 

ii>>)|>p(^^) (9) 



where the coefficient a { i is given by 

^sin[.-V^J 

The relationship between the demodulated signals x,{m) and the parallel source signals s t (m) 

can be expressed in matrix format as follows 

X(m) = AS(m) (H) 

where 

JfW = [x ft W,^ / W,,^,W] r (12) 

S(m) = [j 0 (ifi),... f J l (m),...,^. l (m)] r ( l3 ) 
and the mixing matrix is given by 



(14) 



is 



If there are no frequency offsets, A/ = 0 , a,, = 1 (/ = i) and a,, = 0 (/ * i) . Thus, the matrix A 
the identity matrix. Thus, no ICI occurs. If A/*0, a Ll will not equals zero, which 

results in ICI. The value of a u depends on the frequency error e and on (/-/) mod N , where mod 

is modulus after division, instead of directly on / and i. Hence, the mixing matrix A is a circulant 
matrix, and the inverse matrix of A is also a circulant matrix [8]. Also, all the diagonal elements of A 
equivalent. Usually, in practice the frequency offsets are very small. Hence, eD 1 and 

|a,|D |«,,|,,*/ (15) 

In other words, the absolute values of the diagonal elements in the matrix A are greater than the 



are 



off-diagonal ones. For the wide ranges of the frequency error, -0.4 £ e < 0.4 , and of the sub-carrier 
number, 4 <, N < 52 , the absolute value of the diagonal element and the maximum absolute value of 
the off-diagonal elements are shown in Fig. 2. It is obvious that equ. (15) is generally valid in 
practice. This property is also held for the inverse matrix A, which is shown on Fig. 3. 

3. Proposed ICI/BSS Approach 

Usually, whitening is first applied before the BSS algorithm. X(m) is whitened by a whitening 

matrix Q to produce a white measurement signal vector X(m) , 

X(m) = QX{m) = D*EX{m) (16) 

where D is the diagonal matrix of eigenvalues and E is the matrix of the eigenvectors of the 
covariance matrix # w of X(m) . 

After whitening, the mixing model is given by 

X{m) = QAS{m)=AS(m) (17) 

where the new mixing matrix A is orthonormal. 

To estimate an orthonormal mixing matrix greatly reduces the complexity of the BSS algorithm 
afterwards [8]. Then, the relative gradient method is used to determine the separating matrix W , 
which transforms the white measurement signals i ; (m) into the estimates of the source signals S,(m) . 
The update rule for the separating matrix W is given by [10] 

W = W + ^\l-E^f{s{r))S{r) H ^W (18) 

where n is the convergence factor, I is the identity matrix, E{ ) denotes expectation, and the 
superscript H denotes conjugate transpose. The.estimate source signal vector 

S(m) = WX(m) = [^m),...,l(m),...,s„_ l (m)] T (19) 

and the nonlinear function 

/(S(«)) = [/.(i,(i»)),-./(5,("»)).-./(5».,('»))] r (20) 

The advantage of the relative gradient method is its equivarient property, namely, the 
performance of the relative gradient method does not depend on the mixing matrix [11]. 

After each iteration of the relative gradient method, the separating matrix W is made to be 
orthonormal by the symmetric orthonormalization method [9]. 

W = w{w"w) h (21) 

The separating matrix corresponding to the original matrix A is given by 

W = WQ (22) 

The BSS type algorithm is associated with the permutation and gain ambiguities, which need to 
be resolved in the application. In other words, the separating matrix W is related to the mixing 
matrix A by 

WA = KP 



where P is the permutation matrix and A is the diagnose matrix given by 

(24) 



'A - 0 ' 



A = 



0 ... Av-, 



where gi is the unknown complex number. 

The technique to solve the permutation ambiguity is given as follows. If W is without the 
permutation ambiguity, the separating matrix W is also in same format as the mixing matrix A, 
namely, 

W D h./|.''*' (25) 

The row permutation is applied to the separating matrix W to change its format to satisfy equ. 
(25). In this way, the permutation is solved successfully, and the compensated source signals are 
given by 

S(m) = WX(m) ( 26 ) 

where the parallel compensated signal vector 

S(m) = [s 0 (m),..J l (m),...,s N .,(m)] T (27) 

To resolve the gain ambiguity, additional information is needed to be embedded into the data 
streams. This information is usually in the form of predefined symbol sequences, which is used at 
the receiver detector to solve the gain ambiguity [9,12]. The whole ICI/BSS compensation scheme is 
summarized in Table 1 . 

4. Simulation Results 

The performance of the proposed algorithm is evaluated through the computer simulations. In 
the simulated OFDM system, the number of sub-carriers N equals 8. For each of the parallel source 
signals Sj(m), 30,000 samples of the baseband QPSK modulated signals are used, which are the 
statistically independent. In the relative gradient algorithm, the convergence factor M equals 0.04. 

Also, the nonlinear function /(>>) = \yf y is a good choice, since the typical baseband digital 

communication signals are complex with a negative kurtosis [13]. The effectiveness of ICI/BSS is 
demonstrated by comparing SER of the compensated signal i(n) with SER of the uncompensated 
signal s'(n) , which is obtained from x,(m) by parallel to serial operation. 

For lOdB signal-to-noise ratio (SNR), SER of s(") and s\n) vs the frequency error is given as 
Fig 4. It is obvious the performance degradation due to the frequency offsets from the SER curve of 
s\n) . At the point e = 0 , SER of and s'(n) are equivalent, which means ICI/BSS does not 
degrade the system performance if there are not the frequency offsets. From e = 0 to e = 0.2, SER 
of is almost constant, which means the effect of 1CI is successfully removed after the 

compensation of ICI/BSS. 

Under the fixed frequency offset e = 0.1, SER of s» and s\») vs SNR is given as Fig. 3. 
Thus, under the different noise levels, ICI/BSS is still effective and can greatly improve the system 
performance. 



5. Conclusions 

In this paper the ICI model due to the frequency offsets for the perfect Nyquist channel m 
OFDM is derived. A compensation scheme for ICI in OFDM employing the BSS algorithm, namely, 
ICI/BSS is proposed. A practical technique is presented to solve the permutation ambiguity. This 
approach utilizes the bandwidth efficiently and does not require any prior information about the 
system The simulation results confirmed the excellent properties of ICI/BSS in compensating the 
frequency offsets over a wide range. Also, the performance was maintained under different noise 
levels. 
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Fig.l: Structure of a simplified OFDM communication system 




O The absolute value of the diagonal element 
* The maximum absolute value of the off-diagonal elements 
Fig 2- The absolute value of the diagonal and the maximum absolute value of the off-diagonal 
elements of the mixing matrix A, for the frequency error -0.4 < e < 0.4 , and the sub-carrier number 

4< 7V<52 




O The absolute value of the diagonal element 
* The maximum absolute value of the off-diagonal elements 
Fig 3- The absolute value of the diagonal and the maximum absolute value of the off-diagonal 
elements of the inverse matrix of the mixing matrix A, for the frequency error -0.4 < e < 0.4 , and the 

sub-carrier number 4<;V<52 



Step 1. Whiten of X(m) through a linear transform, equ. (16). 
Step 2. Initialize the matrix W as a random matrix. 
Step 3. Update W by the relative gradient algorithm to, equ.( 18). 
Step 4. Orthonormalize W , equ.(21). 

Step 5. Check the convergence of W , and if the convergence is not reached go back 

to Step 3, otherwise finish the iterations. 

Step 6. Obtain the separating matrix W, equ.(22). 

Step 7. Solve the permutation ambiguity of the separation matrix W. 

Step 8. Obtain the compensated signal vector, equ.(26). 



Table 1. Outline of the proposed ICI/BSS approach 
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Fig.4: SERs vs. the frequency error 
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Fig.5: SERs vs. SNR 



